Spin magnetotransport in two-dimensional hole systems 
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Spin current of two-dimensional holes occupying the ground-state subband in an asymmetric 
quantum well and interacting with static disorder potential is calculated in the presence of a weak 
magnetic field H perpendicular to the well plane. Both spin-orbit coupling and Zeeman coupling are 
taken into account. It is shown that the applied electric field excites both the transverse (spin-Hall) 
and diagonal spin currents, the latter changes its sign at a finite H and becomes greater than the 
spin-Hall current as H increases. The effective spin-Hall conductivity introduced to describe the 
spin response in Hall bars is considerably enhanced by the magnetic field in the case of weak disorder 
and demonstrates a non-monotonic dependence on H . 
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One of the most challenging problems in the physics of 
two-dimensional (2D) electron systems is the excitation 
of spin currents by an electric field E directed in the 2D 
plane. This phenomenon exists owing to the spin-orbit 
(SO) interaction, which brings spin-dependent terms to 
the Hamiltonian of free electrons (intrinsic SO coupling) 
as well as spin-dependent corrections to the scattering 
potential (extrinsic SO coupling); see Ref. 1 for a re- 
view. In the quantum wells grown along [001] crystallo- 
graphic direction in cubic crystals of zinc-blende type, the 
symmetry allows both non-diagonal (perpendicular to E) 
and diagonal (parallel to E) currents of z-polarized spins. 
The diagonal spin current does not appear if the SO split- 
ting is isotropic. The presence of non-diagonal spin cur- 
rent leads to the spin Hall effect, which has been detected 
both in electron 2,3 and hole 4,5 systems by observing spin 
accumulation near the sample boundaries. The original 
theoretical proposal 6 of the intrinsic spin Hall effect has 
been based on the Rashba Hamiltonian 7 describing the 
SO coupling in electron systems due to structural inver- 
sion asymmetry. However, theoretical calculations 8 have 
proved the absence of static intrinsic spin currents for 
this case, and this statement remains true when a mag- 
netic field is applied to the system. 9 Consideration of the 
equation of motion for the spin density operator 10 allows 
one to extend applicability of this result of Refs. 8 and 
9 to any electron system described by the p-linear SO 
coupling. On the other hand, the static intrinsic spin 
currents are present in 2D hole systems described by the 
effective p-cubic SO coupling Hamiltonian, 11 as demon- 
strated by theoretical studies. 12-16 It is believed that the 
experimentally observed spin Hall effect for 2D holes 4 ' 5 
is of the intrinsic origin. 

In spite of the fact that experiments on spin excitation 
are often carried out in the presence of magnetic fields, 
the theoretical work devoted to the influence of a mag- 
netic field on the electric-field-induced spin currents is 
very limited. 17 For 2D hole systems, where the intrinsic 
spin currents exist, the spin-Hall conductivity has been 
calculated 18,19 using Kubo formalism in the collisionlcss 
approximation. Although this approach allows one to 
study the case of strong magnetic field and to describe 
Shubnikov-de Haas oscillations of spin conductivity, 19 



it cannot be used for weak magnetic fields, when the 
cyclotron frequency is comparable to or less than the 
momentum relaxation rate due to scattering. In par- 
ticular, the zero magnetic field limit of the spin-Hall 
conductivity appears to be singular in the collisionless 
approximation. 18 To describe the region of weak mag- 
netic fields, which is the most important for experimen- 
tal studies, a consideration of spin conductivity in the 
presence of scattering is necessary. 

In this Brief Report, the intrinsic spin current is cal- 
culated for the case of a classically weak magnetic field 
H perpendicular to the 2D layer. The elastic scatter- 
ing of carriers is taken into account. The general ap- 
proach to the problem assuming arbitrary SO coupling 
Hamiltonian is followed by application to 2D hole sys- 
tems described by the p-cubic SO coupling. The ana- 
lytical solution obtained below shows that the spin-Hall 
conductivity increases at small H and has a maximum in 
the region where the cyclotron frequency is smaller than 
the relaxation rate. Moreover, it is found that the di- 
agonal component of the spin conductivity appears. For 
this reason, the spin response in Hall bars should be de- 
scribed by the effective spin-Hall conductivity which is 
a combination of diagonal and non-diagonal components 
of the spin conductivity tensor. 

Consider 2D quasiparticles whose states are doubly de- 
generate in spin if the magnetic field and SO interaction 
are absent. In the presence of both the magnetic field 
and SO interaction, the single-particle Hamiltonian is 



H = - \-h„ + V r , 

Am 



(1) 



where tt = ~ih\7 r — (e/c)A r is the operator of momen- 
tum in the 2D plane (xy) and A r is the vector potential 
describing the magnetic field according to H = [V r x A r ]. 
Next, m and e are the effective mass and electric charge 
of the quasiparticles, c is the velocity of light, V r is the 
scattering potential, and <x is the vector of Pauli matri- 
ces. The 2x2 matrix term h n describes both intrinsic 
SO coupling and Zeeman interaction. The extrinsic SO 
coupling effects arc not considered. 

The limit of weak (classical) magnetic field corresponds 
to the condition huj c <C e, where lu c — \e\H/mc is the 
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cyclotron frequency and s is the mean kinetic energy of 
quasiparticles. It is convenient to describe transport phe- 
nomena by using the kinetic equation for the Wigner dis- 
tribution function p prt , which is a 2 x 2 matrix over the 
spin indices (see, for example, Refs. 20 and 21) and de- 
pends on the 2D momentum p, coordinate r, and time t. 
For the spatially-homogeneous and static case considered 
below, the coordinate and time dependence is omitted. 
Searching for the linear response to the applied electric 
field E, one represents the distribution function in the 
form p p = f p eq ^ + / p , where / p is the non-equilibrium 
part satisfying the linearized kinetic equation 20,21 



^■pj /p 



dp 



■+eE-^— = J(/|p), 



(2) 

u(p) = v p + dhp/dp is the group velocity in the pres- 
ence of spin-orbit interaction, v p = p/m, and {a, 6} = 
(ab+ba) /2 denotes the symmetrized matrix product. The 
collision integral J describing the elastic scattering is con- 
sidered under the assumption that the relaxation rate is 
small in comparison to e/h. Assuming that the spin- 
splitting energy 2?i|r2 p | is small in comparison to 2e, it 
is convenient to expand the collision integral in series of 
?3.|S~2 p | /e". Retaining the first-order terms in this expan- 
sion, one obtains 21,22 



J(J\P) 



2tt 



dp' 



(27T?i) 2 ]P - P ' ] 

-{(hp- V)>(/p' 



(/p' - 
-/p)} 



d8(e p i - 



de p > 



,(3) 



where w\ p - p '\ is the spatial Fourier transform of the 
correlation function of the scattering potential V ri and 
£ p = p 2 /2m is the kinetic energy. The condition 
7i|Jlp| <C s also allows one to write the equilibrium dis- 
tribution function in the form f p eq ^ ~ f Sp + h p f' Ep , where 
f' s = df Ep I de p and f e is the Fermi distribution function. 

The calculations below are done in the approximation 
of short-range scattering potential, when u>| p _ p / ~ w is 
constant. It is convenient to use the spin- vector represen- 



tation f p 



f° - 
p 



& ■ f p . In the leading order with respect 



to the small parameter ?i|Jl p |/e p one finds 



f° G f 



r _ ^(E-v p )-^ c [Exvp]-n 

U p _ 6 9,9 ' \ ) 



The components of the vector-function R p are 



d 



^p = ^^-PpGp-o«G p )/;] 



+he[E x n] 



<9e_ 

dp^ f Sp 

d<p p 2 



hu> c 



dip 



JGpf ep 



de r , 
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Equations (5) and (6) are valid for arbitrary f2 p . The 
first and the second terms in Eq. (5) describe spin pre- 
cession and cyclotron motion, respectively. The term R p 
describes excitation of spins by the electric field. Solution 
of Eq. (5) determines the non-equilibrium spin current 
density conventionally defined as 



dp 



[2-Kh) 



r Tr({<T,n 7 (p)}/ p ), 



(7) 



where S = 1/2 for electrons and S = 3/2 for holes in the 
ground-state subband. Applying Eqs. (5)-(7) to electron 
system described by a p-linear SO coupling Hamiltonian, 
one can find that the spin current is zero in the absence 
of Zeeman coupling, as expected. 9 Let us consider the 
2D holes near the bottom of the ground-state subband 
in an asymmetric [001]-grown quantum well. By adding 
the Zeeman term —Suj h <j z to the effective p-cubic SO 
coupling Hamiltonian 11 describing the 2D holes, one gets 



ftp = (up cos 3f, up sin 3ip, 



3w ff /2), (8) 



where k is a constant determined by the Luttinger param- 
eters and confinement potential. Next, u H = gpL B H/%, 
where g is the effective g- factor of holes, and ps = 
\e\h/2moC is the Bohr magneton. Expanding f p in sc- 
ries of angular harmonics, one can solve Eq. (5) exactly 
and obtain the result 



q 7 = (0, 0, q 7 ), q 1 = E^-E/j. 



(9) 



The tensor of spin conductivity, E 7/ 3, contains two con- 
tributions: 



E 7( 3 - E^j + E 7/3 



(i) 



(10) 



/ de p (-f' £p ), (11) 
Jo 



where v = mw/h is the relaxation rate and n = (0, 0, 1) 
is the unit vector in the direction of H. Here and below, 
e is taken to be positive, since the theory will be applied 
to hole systems. With the use of Eq. (4), the equation 



for the vector-function f p is written as 



d 

2[rj p x f p ] - uj c — fp 



Kfp-U (5) 



where (p is the angle of the vector p, and the line over a 
function denotes the angular average (27r) _1 dip . . . . 



- 96 [ dr ( f 1R P [l6 ^ 



-ILU C ) 2 (V + V_+A 2 p ) 

v± = v — iuo c ± 3i(uj c — uoh)- 



,(12) 



In these equations, e zl p is the antisymmetric unit tensor 
and HA p = 2hnp 3 is the SO splitting energy. Notice the 
symmetry properties Y xy = —Y yx and Y xx = E TO . The 

contribution E^ exists owing to equilibrium spin polar- 
ization by the magnetic field in the presence of Zeeman 
coupling. Indeed, the equilibrium spin density is given 
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by s z = S J ( 2 ^)^ r ("' z /p 6g ^) — QuHm/Aizh, and the re- independent of the spin splitting: 

lated spin conductivity is = (s z /erih)cr 7 (3, where <r 7/ 3 
is the Drude conductivity tensor and nu is the hole den- 
sity. In contrast, the contribution E^ is caused by the 
SO coupling. In zero magnetic field, when ui c = luh = 



S (i) 



9e v 2 

%^flV 2 + LU 2 



1 + 



(13) 



0, E^g disappears and E^g contains only non-diagonal 

(Hall) component, E^ = QeT/^fi (see Refs. 12-16), 
where the factor T describes suppression of the spin-Hall 
conductivity by the disorder. In the limit of low tempera- 
ture, when the hole gas is degenerate, T = [l + (i//A) 2 ] _1 , 
where HA = hA PF is the SO splitting energy at the Fermi 
surface and Pf = \J2raep is the Fermi momentum. The 
transition to the low-temperature limit in Eqs. (11) and 
(12) implies — f' e = S(e p — Sf), so the integral in Eq. 
(11) is equal to unity and the integration over energy in 
Eq. (12) is reduced to the substitution A p — ► A. 
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FIG. 1: The magnetic-field dependence of the spin conduc- 

(solid) and of their SO parts 
(dash) for degenerate hole gas. 



tivity components T, xy and T, xx 



and 



The main features of the spin conductivity in the mag- 
netic field are the presence of both non-diagonal and di- 
agonal components and the unusual (non-Drude) depen- 
dence of these components on the ratio u) c /v. Figure 1 
shows the plots of T, xy and E xx as functions of this ra- 
tio for several values of A jv. Notice that the plots for 
A/v = 4 correspond to experimental conditions of Rcf. 
5. The ratio u)h/oj c = gm/2m is estimated as 0.06, us- 
ing g = 0.44 and m/m ~ 0.27 given in Rcf. 19. Since 
this ratio is small, the influence of Zceman coupling on 
T^xy and E^J is weak. However, E^ is considerably sup- 



pressed with the increase of lo c , while E^ saturates at 
the value (9e/87rh)(gm/mo) determined by the g-factor 
and effective mass. Therefore, the contribution E^ dom- 
inates in T, xy at larger oj c / v, especially for the dirty case 
A < v, when the contribution E^ is suppressed. It is 
remarkable that T, xy is a non-monotonic function of the 
magnetic field and increases in the region of weak fields. 
This behavior takes place in the clean regime A 3> v, 
when the spin splitting is not suppressed by the disor- 



der. In the limit A » i/, 3oj c , the contribution E 



(i) 

xy 



According to this expression, the maximum Ely is ap- 
proximately twice larger than the zero-field spin-Hall 
conductivity Qe/SirTi. As follows from the Onsager sym- 
metry principle, T, xy is symmetric with respect to the 
magnetic field reversal. In contrast, the diagonal com- 
ponent Yj xx is antisymmetric in H. This component is 
absent at H = for the particular case of p-cubic SO cou- 
pling considered here. As seen from Fig. 1, E xx changes 
its sign at a finite magnetic field and is strongly sup- 
pressed by the disorder, though the magnetic-field sup- 
pression of this component is weak. The behavior of E 
in the limit A ^> v, 3co c is described by the expression 



(i) 



E (i) 



9e 



VUJ C , 



V 2 + UJ 2 . 



2-3 



luh 6w c (w c - ui H ) 



V 2 + LJ 2 . 



• (14) 




FIG. 2: The magnetic-field dependence of the effective spin- 
Hall conductivity E// for the Hall bar geometry. The dashed 
line corresponds to the limit A»i/, 3lo c given by Eq. (16). 

In the experiments using the Hall bar geometry, the 
electric field acting on the carriers has both longitudinal 
and transverse (Hall) components, the latter is deter- 
mined by the requirement of zero electric current in the 
transverse direction. If the field E y is applied along the 
bar, the Hall field is E x = — (<J xy /a xx )E y = —(io c /v)E y , 
and the spin-Hall current q. x = Y>hE v is described by the 
effective spin-Hall conductivity 



-Y — — V 

- L-i xy Zj a 

v 



9e 



-Rc- 



A 2 v_ 



(15) 



87r?i (y— i(jj c )(y + V- +A 2 ) 

written here for the case of degenerate hole gas. The 
Zeeman part E^ does not contribute to E#. The elec- 



imci 



is 



trie current through the Hall bar is j y = (JqE. 
(To = e 2 nh/mv is the Drude conductivity at H = 0. 
Therefore, the quantity T, H /a determines the ratio of 
spin-Hall current to electric current in Hall bars. Figure 
2 shows a considerable increase of E# with respect to 
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its zero-field value and a non-monotonic behavior in the 
clean regime. If A 3> v, 3cj c , one finds a simple expression 



9e 
8^ 



1 



3w c (w c - u H ) 



(16) 



The Hall bar also carries the longitudinal spin current 
q y = Y,r>E y expressed through the effective diagonal spin 
conductivity S_d = ^xx + {u c /v)Y< xy . 

Since the approach used in this Brief Report does not 
describe the Shubnikov-de Haas oscillations, the results 
presented above can be directly applied if these oscilla- 
tions are suppressed either by the temperature or by the 
disorder. In the general case, Eqs. (10)-(16) should be 
treated as the expressions for the slow envelope part of 
the oscillating spin conductivity. If the relaxation rate 
v is aimed to zero, T^y disappears at a finite magnetic 
field, in agreement with the result of Ref. 18. The effec- 
tive spin-Hall conductivity however, remains finite 
in this limit. 

To conclude, in the presence of a perpendicular mag- 
netic field the spin current of 2D holes is described by the 
spin conductivity tensor containing both diagonal (Yi X x) 
and non-diagonal (E xy ) components. These components 
are non-monotonic functions of H and increase in the 
low-field region. The enhancement of the spin-Hall com- 
ponent Yi X y in comparison to its zero-field value and the 



appearance of the diagonal component T, xx are explained 
by the asymmtery introduced by the Lorentz force, when 
the Fermi surface is shifted in the direction determined 
by the angle arctan(w c /zy) with respect to E; see Eq. 
(4). As the magnetic field increases and oj c exceeds v, 
the components £ xx and Y, xy are suppressed by the field. 
The consideration of the component H xy alone is not suf- 
ficient for description of the spin-Hall response in a Hall 
bar. It is necessary to introduce the effective spin-Hall 
conductivity Sjj which determines the transverse spin 
current q x = Y.hE v proportional to the applied longitu- 
dinal electric field E y . Since the spin density accumu- 
lation near the edges of the Hall bar is estimated as 5 
s e z dae ~ q x m/pp, the behavior of can be directly in- 
vestigated by measuring the magnetic-field dependence 
of this accumulation. The theory suggests (see Fig. 2) 
that Tih is considerably enhanced by the magnetic field 
in the clean systems, where the spin coherence is not sup- 
pressed by the disorder. This condition is attainable in 
the existing samples 4 ' 5 and can be improved by increasing 
the hole density, because A oc ru/ 2 . If A/v ~ 4 (see Ref. 
5), Yih has a maximum at lo c ~ v, which corresponds 
to H ~ 2.8 T if one uses the parameters m — 0.27 tuq 
and hv — 1.2 meV typical for 2D holes in GaAs quan- 
tum wells. 4 ' 5 Therefore, experimental verification of the 
theoretical results is possible and desirable. 
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